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Abstract 



We compute the spectrum of gluonic screening-masses in the ++ channel of quenched 
3d Yang-Mills theory near the phase-transition. Our finite-temperature lattice simu- 
lations are performed at scaling region, using state-of-art techniques for thermaliza- 
tion and spectroscopy, which allows for thorough data extrapolations to thermody- 
' namic limit. Ratios among mass-excitations with the same quantum numbers on the 

gauge theory, 2d Ising and A</> 4 models are compared, resulting in a nice agreement 
with predictions from universality In addition, a gauge-to-scalar mapping, previ- 
ously employed to fit QCD Green's functions at deep IR, is verified to dynamically 
l ■ describe these universal spectroscopic patterns. 
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1 Introduction 



The lattice formulation of gauge theories [T] as QCD allows for ab initio sys- 
tematically well-controlled numerical investigations of most theoretical key as- 
pects, like confinement and chiral symmetry breaking [2J, at non-perturbative 
regime. For instance, at high temperatures or densities when deconfinement 
phase-transition takes place in non-abelian gauge theories, Monte Carlo simu- 
lations have been essential to unveil how matter behaves under extreme ther- 
modynamical conditions [3]. 

Hence, quantum field theories (QFT) on the lattice can be understood as 
classical models of statistical mechanics. Some enlightening interconnections 
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usually appear from finite temperature studies [I], as in the case of a long- 
standing conjecture — designed upon general symmetry arguments — from 
Svetistky and Jaffe [5J relating universal properties from finite-temperature 
SU(N) gauge-theories in (d + 1) dimensions to d— dimensional Z^ — spin sys- 
tems at criticality. This has been corroborated by computations of critical 
(static and dynamic) exponents, universal amplitudes, and correlation func- 
tions: see for example, respectively [6j|7j[8j. 

Still, under the spontaneous breaking of global center-symmetry, SU(N) 
gauge theories were shown — in 4-dimensions and N < 3 — to exhibit a 
rich spectrum of gluonic screening-masses, whose ratios among equivalent ex- 
citations are shared by Z^— spin systems [9lll0| . Despite the effectiveness of 
universality-based reasoning, investigations on symplectic gauge theories com- 
pel for a still elusive dynamical description of deconfinement instead of 
relying completely on Z^— invariance as a physical fingerprint. 

New perspectives have appeared over the last years, such as a unified pic- 
ture of confinement emerging from various analysis of Green's functions of 
non-abelian gauge theories [2"|12|I13|I14|IT5] . While analytical (decoupling) so- 
lutions of Dyson-Schwinger equations [T27T3] and gauge-to-scalar mappings 
|15| predict non- vanishing gluon propagators at deep infrared for d > 2, lat- 
tice computations provide evidences of gluons having dynamically-generated 
masses [14J. 

These features stay unchanged at finite-temperatures, which supports a strongly- 
coupled plasma at the vicinity of phase-transition, thus agreeing with dimensional- 
reduced theories [16] and (universal) results from AdS/QCD correspondence 
[17J. Additionally, properties derived from gauge-dependent gluon propaga- 
tors would also be compared with spectroscopical analysis employing gauge- 
invariant correlators [TS] (free from Gribov ambiguities [19]), which may deepen 
physical insights into QFT- thermodynamics. For a more comprehensive re- 
view on gluonic bound-states using methods from spectroscopy, at zero and 
finite-temperatures, see also |20| . 

In this article, we investigate through Monte Carlo simulations the gluonic 
screening-spectrum of quenched three-dimensional Yang-Mills (YM) theory on 
the lattice. State-of-art techniques for spectroscopy, such as the variational- 
method [21J and recursively-smeared interpolators |9]fT0] . are employed to com- 
pute excitations on the ++ channel. Using an improved thermalization algo- 
rithm [22] we carry out simulations at asymptotic region with ameliorated 
critical-slowing-down (CSD) |23| . UV-discretization and IR-cutoff effects |24| . 
We compare mass-ratios of YM theory and exactly-solvable bidimensional sys- 
tems in the same universality- class (i.e. the Ising and A0 4 models) [5j, therefore 
enabling us to check predictions from universality [5,9, 10J and other dynamical 
mechanisms |11|IT5] . 
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This work is organized as follows: in section II there is a review of the major 
aspects of universality and its implications for thermal excitations on gauge 
theories and spin-systems at the vicinity of critical points. Section III outlines 
the general simulation setup, as well as the modified heat-bath algorithm 
used for ensuring more efficient thermalization of gauge fields. Spectroscopical 
techniques for extracting screening masses are the topic of Section IV. Section 
V presents numerical results and comparative data analysis. Conclusions of 
the present work and its outlook are discussed in Section VI. 



2 Universality: field theories and spin-systems 

The finite-temperature formalism for d— dimensional gauge theories on the 
lattice uses an asymmetric euclidean spacetime, whose physical volume 
V = a d ■ Lff 1 • Lt has L$ and Lt as spatial and temporal (adimensional) lat- 
tice sides, respectively. Thus, under the assumption a ^ (a ■ Lt) < {a ■ Lg), 
where a is the (dimensionful) lattice spacing, one can set the temperature of 
physical equilibrium as T _1 = a ■ Lt- As a consequence of the formalism, in 
the vicinity of phase transitions, arguments from universality [2] state that 
finite-temperature QFTs for SU(N) gauge groups in (d + 1) dimensions be- 
long to the same universality class of globally ^TV-invariant spin-systems in d 
dimensions. 

For the particular case of SU(N = 2) gauge theory, a second-order (i.e. crit- 
ical) deconfinement phase-transition is expected. This theory can be imple- 
mented on the lattice by the Wilson action 

Sw [U] = /?£ £ {l - \Re{TrP, v )\ , (1) 

using SU(2) gauge-links (x) and the plaquette 

= (x) U u (x + fia) U' 1 (x + va) V~ x (x) , (2) 

with lattice-coupling = 2N / g 2 s a A ~ d and gauge-field coupling constant g 2 s . The 
action Eq. ([T|) is invariant under transformations generated by center-elements 
of the gauge-group, hence, spatial loops — in contradistinction to temporal 
ones — are insensitive to temperature-induced symmetry breaking. Then, a 
natural order parameter is given by a temporal loop, the so called Polyakov 
loop L = (L (x,y, z)) = (T r n™Zi T Ut (x,y, z,an)*j [25], which measures the 
potential energy necessary to free a quark [3]. 

On the other hand, the 3d Ising model is a well-known spin-system that un- 
dergoes a second-order phase-transition, for (3i = (5 'critical, due to a global 
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Z 2 — symmetry breaking universally related to the finite-temperature (3 + l)d 
YM theory. It is described by the action 

Si = -Pi s n s m- (3) 

<n,m> 

In this same universality class is the 3d lattice-regularized (p 4 theory, whose 
action is 



S* = -P* E 0n0 m + E^ + AE(^- 1 ) 2 ' ( 4 ) 

<n,m> n n 

In both cases, the order parameters are expectation values of fundamental 
fields, e.g. the magnetization (m) or the v.e.v. ((fj of the scalar-field, which 

behaves as the Polyakov loop (^Lj in the YM theory [9]. 

The analogy among order-parameters m, (f) and L can also be extended to 
their correlation functions. The connected correlation function Gq is usually 

expressed as 

G (|r 2 - n|) c =< O (r 2 ) Ot (n) > - < O (r 2 ) >< Ot (n) >= 

which is a sum over each mass (rrii) of the spectrum of the operator O, in an 
L$— periodic lattice. Thus, when O = m, one gets excitations from a magnetic 
system or, with O = L, the spectrum of screening-masses for the gluonic field 
is recovered. While in the first case the spectrum is generated by a magnetic 
phase-transition that breaks the global Z 2 — symmetry, in the latter case it is 
due to a spontaneous .^-center symmetry breaking of the SU(2) gauge group 
[5|9|10j . 

Although inherently different in nature, the aforementioned phase-transitions 
constitute critical phenomena in the same universality class [4J. Thus, one 
may expect that universality [5] can predict some dynamical aspects of gauge 
theories, such as mass-ratios |2S] of excitations in Eq.flS]), to be shared with 
statistical-mechanical systems. While evidences supporting this hypothesis are 
accumulating |9~fT0] . arguments entirely based on universality are not enough 
to describe deconfinement in sympletic (or exceptional) gauge theories jX 1 j . 
Thus, a broader dynamical picture may be needed |15] . 

An improved understanding of such dynamics may be obtained by compar- 
ing the spectrum of gauge theories and exactly- solvable systems in the same 
universality class. For instance, the 2d Ising model is well-described by a free- 
fermion field theory near criticality, whose spectrum of excitations — for n 
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integer — turns to be 

M n = nM , (6) 

where M is a mass-gap proportional to |t|: the deviation from critical tem- 
perature [27]. At the same time, the spectrum of A0 4 theory is analytically cal- 
culated by a gradient-expansion (28], which results with mass-scales {^-broken) 
— and integer k — for the broken phase 

M k = kA broken . (7) 

Therefore, one can apply Monte Carlo methods to exploit whether the screen- 
ing spectrum of the YM theory matches expectations from Eq.(j6]) and Eq.(J7J 
near criticality, a regime where correspondences are expected, but other non- 
perturbative techniques are less effective (since the theory is dynamically triv- 
ial just in 2d; i.e., at the T — > oo limit [29.30J). 



3 Algorithms 

Since in the vicinities of phase-transitions the CSD phenomena [23] afflicts 
the generation of statistically independent gauge-configurations more severely, 
continuous efforts toward the improvement of thermalization algorithms are 
ubiquitous for Monte Carlo simulations. To ensure an efficient thermalization 
of gluonic fields, we applied our modified heat bath algorithm (MHB), which 
was shown to be faster [22J than the usual heat bath (HB) algorithm [3"Tf3"2"] . 

Usually, when an SU(2) lattice-gauge theory is considered |22] , one can fac- 
torize the action Eq.([T]) as a sum over many single-link actions Si_n n k as 

S^ li nk = ~Tr[U fl (x)H fl (x)}, (8) 

with Ufj, (x) G SU (2) ; (x) = A M (x) (x) is a sum over staples — i.e., it 
is proportional to an SU(2) matrix — where (x) G SU (2) and iV M (x) = 

^detH^(x). 

Then, by using Eq.Q and the invariance of the group measure under group 
multiplication [31], the usual HB update is obtained 

U° ld (x) — ► U; ew (x) = VHl (x) , (9) 

where V = VqI + i-v ■ a G SU (2) is randomly generated by choosing v o 
distributed as y/1 — v 2 exp (flNvo) cLvq and v is randomly chosen in M. 3 . Our 
algorithm |22| proceeds as the HB algorithm to generate the updating matrix 
V, except for the additional step 
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• Transform the new vector-components of V as v — >■ —sgn(v ■ w)v, 

where W = w I + i ■ w ■ 3 = U° ld (x) (x), and sgn is the sign function. 

Still, MHB may be seen as a modification of the overheat-bath algorithm (OH) 
[33j that incorporates a micro-canonical move [34j into a heat-bath step. How- 
ever, while in MHB all components of the vector v are randomly generated 
(i.e., except for its sign) in very close similarity to HB, the OH algorithm 

deterministically sets v = —w (and renormalizes it as \\v\\ = \Jl — Vq). Thus, 
OH incorporates a microcanonical move in an exact, but maybe non-ergodic 
(see for example discussions in [22j) algorithm, which is corrected — by con- 
struction — in the MHB version. 



4 Spectroscopical methods 



The spectrum of excitations of a field theory on the lattice can be directly 
extracted by brute- force least-squares fitting from Eq.([5]) using a multiple- 
exponential decay function. Despite being straightforward, this method has 
several drawbacks and allows for reliable results just when high-quality statis- 
tics is available. More robust techniques such as Bayesian fit [35J or the Evo- 
lutionary fitting |36| would constitute alternatives to overcoming these diffi- 
culties. However, here we use another state-of-art approach, namely the vari- 
ational method |21| . broadly employed in hadron spectroscopy |37j . 

On the variational method a proper set of base-operators (i.e., the interpola- 
tors Oi) has to be chosen for building a cross-correlation matrix 

dj (\r\) =< % (r) O) (0) > - < Oi (r) >< O) (0) >, (10) 

which may be diagonalized in a generalized eigenvalue problem (properly nor- 
malized at some r < Ls), resulting in 



C (r) • v n = X n (r, r ) • C (r ) • v r , 
where eigenvalues behave as 



A„ (r, r ) a e 



-(f-foK 



l + 0(e 



-(r-r )Am r 



(11) 



(12) 



Generally, Am n is the mass difference to the closest lying state, where each 
interpolator Oi (projected to defined momentum-states) has quantum numbers 
in a given channel. 

Within a large enough basis of independent interpolators, each eigenvalue of 
Eq. (fTU|) will decay as the leading order of Eq.( lT2]) . Thus, the slowest decay- 
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mode of eigenvalues is associated with the ground state, while the fastest 
one gives the highest-excitated state. This implies that simple fits of a few 
parameters may be nicely done, since single-exponentials dominate the signal 
in the whole range. In order to identify proper fit-ranges, it is useful to locate 
stable plateaus not only for effective-masses of X n 



f„" (r + 1/2) = In (- 



<"< r+1 '^ h ImTi)' (13) 



but also in their associated eigenvectors v~k, which work as fingerprints of each 
state. 

We are interested in a set of interpolators that generates scalar screening- 
excitations in the ++ channel. An efficient method for building this is to 
apply recursive smearing-steps over usual wall-averaged (i.e. zero-momentum) 
Polyakov loop operator |9], which is defined by 

L(x) = -— £ L(x,na). (14) 

\ L S)y n= l 

It enables access to all length-scales on the lattice through the computation of 
cross-correlation functions among nth-smeared interpolators (i.e. p( n >), which 
are constructed iteractively (for increasing smearing steps n) from the usual 
Polyakov loop (L (x)) . That operators are assembled following the rule 

p(°) ( x ) = L(x) 

p(n+l) ( x j = gign ^ bj _ w ) | w | + u ^p(n)^j (15) 

u = \ fpW (y - a) + P( n) (y + a) 
where u G (0, 1] and we have taken u = 0.1. 

In addition, a compromise between signal-to-noise and maximal linear inde- 
pendence of correlators must be attained |38j, which is achieved by looking 
for a set of interpolators that minimizes the conditioning number (k) of the 
normalized correlation matrix (CV,) on slice r = r t 

4- = i aAn) • (16) 

y[c* (r t ) Cjj {r t ) 

Thus, noisy interpolators are eliminated by considering the signal-to-noise 
ratio of the diagonal elements of Eq. (Tl6|) . Additionally, their independence is 
estimated remembering that n = 1 for a completely orthogonal basis (while 
k — y oo for increasing levels of linear-dependence). 
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5 Numerical results 

In our finite-temperature simulations we have used the Wilson action Eq.flT]) 
and asymmetric lattices L| x Lt, where we take = 8 and L5 = {50, 70, 90}. 
The critical coupling we adopted is known with high accuracy to be /3 = 12.63 
[39f4"5] for this case. 

In 2+1 dimensions the gauge coupling g 2 s has the dimension of mass and sets 
a scale for the theory. For instance, the lattice spacing (a) can be calculated 
using the string-tension (a) — as a numeric physical input — on known /3- 
string tension relations for SU (N) in three dimensions |40j . We have assumed 
yfo = 0.440 GeV as in [4"T|4"7] - i.e. a zero-temperature 4d numerical value 
— which gives a ~ 0.0503 fm for SU(2). 

So, our simulations are in the asymptotic region, where discretization effects 
(UV) due to coarse-lattices [24J would be weak. We employed the MHB ther- 
malization algorithm, which incorporates an overrelaxation step, to generate 
respectively {5M, 10M, 15M} gauge-configurations at critical temperature for 
lattices Ls = {50, 70, 90}. The statistical independence of gauge configurations 
was checked by computing the auto-correlation time r, ni 



for the Polyakov loop operator (i.e. O = L) among successive samples. An au- 
tomatic windowing procedure with c = [4, 10] and the Madras-Sokal formula 
for error estimation were employed [23]. We kept lk independent gauge con- 
figurations for spectroscopical analysis. Statistical error-bars were computed 
by the bootstrap method [T|23] with 70 A; repetitions, see Figure (1). 

The interpolators in Eq. (TT5l) were pruned as a function of their smearing level 
using the aforementioned methods. The better suited iteraction-levels we could 
determine were n = {3, 9, 15, 21, 27} , with k ~ 1.2. Although up to five inter- 
polators were used to build cross-correlation matrices, as in Eq. dlOp . no more 
than three mass-excitations could be recovered by our variational approach. 

To proceed the determination of the spectrum of screening masses, we per- 
formed least-squares global-fits to exponential decays of eigenvalues of cross- 
correlation matrices assuming 



as shown in Figure (1) and Table (1). Near phase transitions, some finite-size 
(and tunneling) effects were reported [9"]fTU] to induce sub-leading contributions 




(e>iO i+k )-{e>i) 2 

(Oi 2 )-(Oi) 2 



EtLiPo(k), 



(17) 




(18) 
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Figure 1. Eigenvalues of the cross-correlation matrix for p( n )-operators with 
n = {3, 9, 15, 21, 27} . Panel A: L 5 = 50, panel B: L s = 70, panel C: L s = 90, 



Mass 


£50 


fit range 


L70 


fit range 


L90 


fit range 


M 3 


0.191(8) 


[3,47] 


0.134(3) 


[2,68] 


0.098(4) 


[3,87] 


M 2 


0.107(3) 


[2,48] 


0.081(2) 


[4,66] 


0.058(2) 


[5,85] 


Mi 


0.065(1) 


[2,48] 


0.0467(5) 


[3,67] 


0.0329(5) 


[3,87] 



Table 1 

The screening-masses in lattice units (Mass) obtained by fitting the eigenvalues of 
the cross-correlation matrix of ^"-correlators by Eq. (ll8p . with n = {3, 9, 15, 21, 27} 
for lattices Ls = {50,70,90}. The ranges shown correspond to best-fits with good- 
ness ( X 2 /dof) L50 = 0.306, (i? 2 ) L50 = 0.993; ( X 2 /dof) L7Q = 0.139, (R 2 ) L70 = 0.997 
and {x 2 /dof) L90 = 0.248, (R 2 ) L70 = 0.993. Error-bars were computed by 70/c 
bootstrap-repetitions performed over 1000 independent configurations. Integrated 
auto-correlation times (ri n t-p) associated to Polyakov loops (i.e. P^) were used 
for checking the statistical independence of samples, they were computed by au- 
tomatic windowing procedure with c = [4, 10] . The error estimation for each 
Ti n t-p is from Madras-Sokal formula [23], resulting in (Tint-p)^ = 5366(645), 
(n nt -p) L70 = 9526(1013) and (r mt _ P ) L90 = 15108 (1402). 



to correlators Eq. (1121) . which are handled by the constant A. Each effective- 
mass is obtained from the best fit inside of plateaus, where (at least three) 
neighbour masses agree within one error-bar. 

We also performed a finite-size scaling (FSS) study using the measured masses 
(m n ) for extraction of their infinite- volume limit [42] . A fitting ansatz m n = 
a n /L + b n suggests b n = for all masses, while m n = a n L~ Un corroborates 
with v n = v, within one standard-deviation. Thus, the following functional 
dependence 

m n = a n L~ u , (19) 
was employed for fitting simultaneously all data points, as is seen in Figure 
(2). 

The obtained exponent v = 1.16(6) shows that at IR limit all mass-excitations 
(m n ) scale as the inverse correlation-length (£), as in the critical 2d Ising-model 
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|26| . The average ratios computed from face values in Table (1), using error- 
propagation formulas, are compatible with results calculated from coefficients 
a n = [ ai = 6.4(1.7), a 2 = 10.7(2.9), a 3 = 18.4(5.0)] from the best-fit (i.e. 
x 2 /dof = 0.808, R 2 = 0.993) with Eq.(lT91. Thus, it is pointed out that 
ratios for screening-masses in ++ channel of the critical finite-temperature 
YM theory in (2 + 1) d are m^/wii = 1.69(15) and m^jm\ = 2.89(35). 

Hence, the measured spectrum can be fitted by b n = 0, a behavior associated 
with a tower of massless excitations is seen at thermodynamic limit. Morever, 
our data agree with universality arguments Eq.([6]) and, so far, the prediction 
Eq.© from mapping [15J stays valid for d ^ 4. 



io'- • • 



50 60 70 80 90 
L 

Figure 2. Scaling of effective-masses as a function of lattice sides (L). Dot-line is the 
best global-fit from Ansatz in Eq.flE]), its goodness is {x 2 /dof = 0.808, R 2 = 0.993). 
The exponent obtained v = 1.16(6) and the (a n ) coefficients imply ++ screening- 
mass ratios compatible with m2/mi = 1.69(15) and m^jm\ = 2.89(35). 



6 Concluding remarks 

We have computed the gluonic spectrum of screening masses (in ++ channel) 
of quenched 3d Yang-Mills theory at critical temperature. State-of-art tech- 
niques from hadron-spectroscopy, such as the variational-method and recursively- 
smeared interpolators, were employed with our new thermalization algorithm, 
to ensure accurate results at the critical region. While the measured spectrum 
is expected to present no significant systematic effects from UV-cutoff — since 
we employed very fine lattices — a noticeable dependence on the spatial lattice 
side (i.e. on the dimensionless quantity Ls) was observed. 

Therefore, in the way to obtain the mass-excitations in the thermodynamic 
limit a FSS extrapolation was performed, which unveiled a tower of mass- 
less excitations, as expected by dynamical |15.28] and universality-based ar- 
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guments [5|9|10f27| (i.e. similar to the spectrum of scalar/free-fermion field 
theories). A result that strongly resembles patterns from bound-states of glue- 
balls analytically calculated for SU(N) theories, at the N — > oo limit |43| . 

Hence, traditional criteria for confinement are based on the gauge-dependent 
behavior of gluon-propagators at the infrared limit [2j ; connecting finite-temperature 
studies in different gauges [16] to gauge-independent [18J results is worthwhile. 
An analytical development in this direction was introduced in [T3], where a 
mapping from SU(N) gauge theories to N— scalar fields consistently fits 4d 
gluon-propagators at deep infrared. Moreover, it agrees with modern findings 
concerning YM theories for more general gauge groups [11], where a subtle dy- 
namical balance among group generators and the size of group center seems 
to lead to occasional universal behavior. 

Whether such assertions are valid for the finite-temperature 3d YM the- 
ory near criticality is a matter for numerical verification. So, we have com- 
puted mass-ratios of SU(2) gauge theory and compared them to ratios from 
universality-related models, namely the A0 4 and Ising-field theories. Our re- 
sults match both universality and dynamical-mapping hypothesis [5f9"|15|l28| 
in Eq.([6]) and Eq.flTJ); which harmonizes with confinement-scenarios [29, 30, 44 J 
presenting (confined) massive gluonic-excitations in d > 2. 

For future research we leave the computation of screening-masses for general- 
ized YM theories [llj using smeared interpolators like Eq. (115|) . though built 
upon order parameters, as the dressed Polyakov-loop, associated with chiral 
Dirac operators |37||4"6] . The resulting spectral patterns arising from such in- 
vestigations could be straightforwardly related to the behavior of usual YM 
propagators by methods presented in [47], which may bring further insights 
into dynamical aspects connecting universality to chiral and deconfinement 
phase-transitions |4"8] . 
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